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Abstract 

Fully packed trails on the square lattice are known to be described, in the long distance limit, by a collection of 
free non compact bosons and symplectic fermions, and thus exhibit some properties reminiscent of Brownian motion, 
like vanishing fuseau exponents. We investigate in this paper the situation for their non-intersection exponents. 
Our approach is purely numerical, and based both on transfer matrix and Monte Carlo calculations. We find some 
evidence for non-intersection exponents given by CFT formulas similar to the Brownian case, albeit slightly different 
in their details. 

1 Introduction 

The study of conformally invariant random curves embedded in two-dimensional space has a long history in theoretical 
physics and probability theory. Recently the subject has gained new momentum due to the advent of a set of rigorous 
methods known under the name of Stochastic Loewner Evolution (SLE). 

Two classes of random curves have been especially well studied. In the first class one finds self-avoiding and mutually 
avoiding curves, which come in several variants depending on their precise microscopic definition, and in particular 
whether they occupy a finite fraction of space in the continuum limit. These curves admit a height representation 
(roughly speaking, using the curves as level lines of the height) and as such are amenable to the use of Coulomb Gas 
(CG) techniques. Also, the corresponding (twisted) vertex models can be studied as quantum spin chains, i.e., using 
Bethe Ansatz (BA). As a result, the bulk critical behaviour of self-avoiding curves is very well understood, whereas 
important questions regarding their surface behaviour still remain open [T] . Many of the exact results obtained by the 
CG route have recently been confirmed, or extended, by the SLE approach. 

The second example of random curves is that of Brownian motion, which is best studied directly in the continuum 
limit. Profound links exist between this and the foregoing case, via the use of SLE [2] or quantum gravity [3] methods. 
To give but one example, the external hull of the Brownian curve turns out to be in the universality class of the 
(standard, dilute) self-avoiding walk J3J. On a microscopic level, however, there are two crucial differences. First, 
a Brownian motion can intersect itself, making the application of height mappings and CG techniques break down. 
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Figure 1: The three allowed vertices in the Brauer loop model, with weights 1, l,x. 




Second, it can pass through the same site an infinite number of times, making it more difficult to use the BA approach^. 

As a result, up to this day, the non-intersection exponents for Brownian motion, which were identified by Duplanticr 
and Kwon [5j [3] twenty years ago, and whose values were rigorously established recently [2], cannot be derived using 
the BA or variants of the CG method applied directly to a suitable discrete lattice model in flat space. 

The present work started with the recent discovery of lattice models that bear a close resemblance to Brownian 
motion, yet are tractable by the usual techniques. These models are what is sometimes called dense trails, and can be 
obtained by allowing self-intersections as well as mutual intersections in a dense loop gas, while still requiring that the 
curves pass through each lattice link (resp. site) once (resp. twice) at the most. The introduction of such intersections 
is well known not to change the universality class in the dilute case, but — surprisingly maybe — does affect it profoundly 
in the dense case (for an early discussion see [5]). It was in particular argued in [7J that for such a model with a loop 
fugacity n < 2, all the fuseau exponents vanish exactly, and the corresponding correlation functions are described by 
non compact bosonic fields, exactly like in the pure Brownian case. This was checked analytically and numerically in 
the case of fully packed trails (also called Brauer loop model) [8j [7] . 

It is of course tempting to ask how close fully packed trails and Brownian motion actually are, and a natural route 
to investigate this question is to study the non-intersection exponents. 

Fully packed trails, or Brauer loop models, depend among other things on the fugacity n of loops. We will 
concentrate here on the case n = (hence, of a single trail), but notice that the case n = 1 has recently been studied 
numerically by Kager and Nienhuis [9]. These authors found that the hull distribution was compatible with that of 
a. reflected Brownian motion. In our work, we define a scale-invariant "escape path" T within a Brauer loop model 
with fugacity n — 0. Our main findings arc that while the fractal dimension of T agrees (marginally) with that of the 
self-avoiding walk, the non- intersection exponents arc definitely different from those of [HI H] , though might well be 
given exactly by a related formula. 

The paper is organized as follows. In section [2] we define precisely the model of fully packed trails to be studied, 
and in section [3] we introduce the escape path T and discuss how it can be made scale invariant. A qualitative phase 
diagram is established in section[4] Finally, the non- intersection exponents and the scaling properties of T are measured 
using transfer matrix (section [5]) and Monte-Carlo (section [S]) techniques. We give our conclusions in section [JJ 

2 Brauer loop model on the square lattice 

We consider, on the square lattice, the Brauer loop model defined by the three vertices a, b and c, shown in figure [1] 
The local Boltzmann weights are given by : 



We give an additional weight n for each closed loop. 

1 We note that there has been some important progress recently [4] in the study of Bethe Ansatz for statistical mechanics systems with 
an infinite number of degrees of freedom per site in relation with quantum spin chains on non compact groups. 
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Figure 2: Rcnormalization flow of the loop model, for the parameter x. 



The renormalization flow for the parameter x is shown in figure [2] The fixed point x = corresponds, for 
— 2 < n < 2, to the low-temperature phase of the 0(n) model [10]. The fixed point a;* attracts, for integer [7] 
and presumably also real values of n, all the points in the phase x > 0. This phase contains an integrable point 
x\nt = (1 — n/2)/2 [8j. At this particular value of x, the i?-matrix (expressed in the connectivity basis) satisfies the 
Yang-Baxter equations for any real value of n. 

The goal of this paper is to study the non-intersection properties of the Brauer loop model. Recall that in the 
Brownian case, the simplest such property is the disconnection exponent, which governs the probability that the origin 
of a single Brownian path remains accessible from infinity without the path being crossed. Let us recall some of the 
well established results in this case. If one considers a random walk of length S, the origin stays connected to infinity 
with probability P(S) ~ S^ 1 * 1 , with d = 1/8. One can also imagine growing the walk until it reaches the circle of 
radius R, and ask the probability with which it has left the origin connected to infinity, in which case P(R) ~ i? -2 ** 1 . 
Finally, one can consider a two-point correlation function defined as the weighted sum [5] over all walks going from 
the origin to a point r for which cither extremity remains connected to infinity. When the monomer fugacity equals 
the critical value, G(r) ~ |r|~ 2 ' Tl with x\ = 2(i = 1/4. 

By analogy, we define an interface as a line (of arbitrary length) that the Brauer loops are not allowed to cross. 
In particular, we adress the following problems : 

1. Define the above loop model on a finite lattice of iV sites. Let the origin be a point inside the lattice, and r a 
point on the boundary of the lattice. Let P(r) be the probability that there exists an interface going from the 
origin to the point r. How does P(r) decay when r = |r| goes to infinity ? 

2. When this interface exists, it is, by definition, a self-avoiding walk — for easy reference we henceforth refer to it 
as r. The original loop system around T acts like an environment, affecting its geometrical properties. What is 
the universality class of T ? 

3 Definition of a fluctuating interface 

It turns out that in the Brauer loop model with a single loop, for any x > 0, the probability P(r) behaves as : 



where £ is a length-scale depending on x. This behaviour is profundly different from the one of Brownian motion, where 
the probability -P(r) for a long Brownian trajectory decays algebraically (a consequence of subadditivity arguments 
|llj). This indicates that, without modification, the interfaces in the Brauer model at n = do not have interesting 
scaling properties. 

Let us see more explicitely why such behaviour arises. Define three partition functions Zq(N), Z'(N\0,r) and 
Z (N\T) as follows. 

• Let Z (N) be the partition sum of the loop model for N sites. 



P(r) ~ exp(-r/0 



(2) 
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• Let Z'(N\Q,r) be the partition sum of the loop model with the constraint that there exists an interface not 
crossed by the loop, going from to r. When a loop configuration admits an interface T, then this interface is 
unique (or else there would be several loops). Thus, one can rewrite the sum Z'(N\0,r) by grouping the terms 
corresponding to the same interface T : 

Z'(N\0,r)= J2 Z ^ N \ T ) ( 3 ) 
resAW(o,r) 

where SAW(0, r) is the set of all self-avoiding walks with ends fixed to and r, and Zq(N\T) is the partition 
sum of the loop model with the interface defined by the self-avoiding walk configuration T. 

Consider the presence of an effective interface energy in Zq(N\T) : 

Z (N\T) ~ Z (N) expH x F a (x)] (4) 

where I is the length of T, and F s is a free energy per unit of length of T. The number of self- avoiding walk configurations 
of length I is approximately (^saw)', where /x$aw = 2.6381(5) [12] is the connectivity constant of self-avoiding walks 
on the square lattice. If the quantity (F s — log/^sAw) is positive, then the behaviour of Z'(N\0,r) is governed by 
configurations where the interface has minimal length : 

I ~ r (5) 
Z'(N\0,t) ~ Z (N) cxp[-(F s -lo gA/S Aw)r] (6) 

This accounts for the exponential behaviour of P(r) = Z'(N\0,r)/Z o (N). 

In order to define a scale-invariant interface, we need to allow its length I to develop big fluctuations. A simple 
way to do this is to introduce an additional Boltzmann weight to , where I is the length of I\ Thus, we are interested 
in the study of the partition function : 

Z(N\0,r)= ]T w l Z (N\T) (7) 
reSAW(o,r) 

This function contains two external parameters : x, the fugacity of a crossing in the loop model defined by Zo(iV|r); 
and w, that controls the length I of the interface T. A typical configuration contributing to the partition sum ([7]) is 
shown in figure [3] 



4 Phase diagram 

Using renormalization group arguments and anticipating on our numerical and analytical results, we try to give a 
consistent picture of the phase diagram of the model defined by equation ([7]) . 

Let us fix the value of x > 0, and observe the system as we vary w. In this picture, the interface Y is a self-avoiding 
walk, with its self-interactions induced by the surrounding loop model. According to equation (T3|) , we expect a phase 
transition at w*(x) = exp[i 7 ' s (a;)]//xsAW- 

Note that it may be possible to check the previous relation. We can compute numerically the value of to*. If we 
could also compute independently the value of F s , then wc would be able to compare these values with the connectivity 
constant fJ^sAW, which is a well-known constant [12j . 

Given w* , the qualitative behaviour of the system as a function of to is expected to be as follows : 
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Figure 3: A typical configuration contributing to the partition sum Z(N\0, r). 

• If w < w*, the behaviour of Z(N\0, r) is governed by configurations where the interface has minimal length : 
this phase is "non-critical" from the point of view of the interface. The statement ([2]) amounts to the inequality : 
w* > 1, for every x > 0. In the phase w < w* , the interface acts as a non-fluctuating boundary for the loop 
model. So we expect the points of this phase to flow towards one of the fixed points described in section [SJ 

• If w = w* , the interface length I fluctuates above its minimal value r, with fluctuations of the order of the system 
size. It is a fractal object, with Hausdorff dimension 1 < df < 2. We call this point the 'T-dilute" critical point. 

• If w > w* , two sccnarii, illustrated in figure[31 arc possible a priori. In scenario 1, there exists a fixed point wd 
such that w* < wd < oo, and the whole phase flows towards it. In scenario 2, no such fixed point exists, and 
the whole phase flows towards the point w = oo. Numerical studies suggest that scenario 1 is correct. 

This behaviour with respect to w is very similar to the one of ordinary self-avoiding walks when driven into the 
dense phase by a large monomer fugacity |13j . The interface T plays the role of the single chain, and the parameter 
w is the monomer fugacity for the chain T. 

To summarize, the phase diagram is given in figure [5] The points (1) and (2) are the critical points of the 
intersecting loop model with a fixed boundary. The point (3) (resp. (4)) is the 'T-dilute" point for x > (resp. 
x = 0). The point (5) (resp. (6)) is the 'T-dense" point for x > (resp. x — 0). The line passing through (3) and (4) 
has equation : w = w*(x). 

Let us focus on the critical point (4). Since x = 0, only vertices of type a and b are allowed, and the point 
(4) is equivalent to the Q-state Potts model on the square lattice, in a particular limit : in the Fortuin-Kasteleyn 
formulation, a cluster has fugacity Q — ► and a bond has fugacity v = \[Q. The leading term of the partition sum 
describes equiprobable spanning trees. In this context, the interface T is exactly the set of red bonds of the spanning 
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Figure 4: Two scenarii for the flow diagram of the parameter w. Numerical studies suggest that scenario 1 is correct. 

tree (the red bonds are the bonds of the spanning tree which disconnect from r if removed). The fractal dimension 
of this object is 5/4 HUE)]. 

The point (2) is just equivalent to the Q-state Potts model in the limit described above, with a fixed cut line 
uncrossed by the loop. 

5 Transfer matrix formulation 

5.1 Mapping from the plane onto the cylinder 

In order to estimate partition sums such as Zq(N), Z(N\0,r) on a cylinder of fixed width L and length M — > oo, 
we compute the leading eigenvalues of the transfer matrix. Denote them A^\ A^\ , ... in decreasing order of 
modulus. Then the eigenvalue defines the free energy density by surface unit in the following way : 

Z%{M,L) ~ (A«) M (8) 
/f^bg^I) ~ (9) 

For conformally invariant models, finite-size corrections to give access to the central charge c and to the operator 
dimensions X® [ME] : 



fL ] = foo + if 2 +o{L-*) (10) 



6i 2 

f?-tf - 2 -^ 1 + o(L^) (11) 
where the dimensions are defined in plane geometry by 



S«(0,r)/^~l/r» W (12) 
5.2 Transfer matrix for the loop model 

Our convention is a transfer matrix Tl acting in the vertical direction, from bottom to top. Let the width L of the 
cylinder be an even integer. Consider the strands living on the vertical edges of the lattice. A row configuration is a 
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Figure 5: Phase diagram of the model with n — 0, in the x — w plane. 



pairing of these strands. There are (L — 1)!! = lx3x5x...(L — 1) pairings of L objects. If a and /3 are two row 
configurations, then the matrix element (Tl)^ is the total Boltzmann weight of L- vertex configurations that map a 
to (3 (see figure |6]). This Boltzmann weight includes a factor n each time a loop is closed. 
On a cylinder of length M, the partition sum is given by : 

Z , cyl (M,L) = (a out | (T L ) M lam) ~ (a£ 0) ) M (13) 

where \cxi n ) and |a ou t) are two particular vectors coding the closing of the loop at each extremity of the cylinder. For 
example, one may choose to close the loop as illustrated in figure [7] 

5.3 Definition of an interface on the cylinder 

Let r be an interface connecting a point A on the bottom circle of the cylinder to a point B on the top circle of the 
cylinder (see figure |5J). Let Zo }Cy \(M, L\T) be the partition sum of the loop model constrained by the interface T. 
Define, as in equation ([7p. the partition sum Z cy \{M 1 L) : 

Z cyl (M, L) = J2 w l Z 0iCyl (M, L\T) (14) 
r 

Note that the vector (Tl) M |a; n ) may contain components on configurations which admit several interfaces. This is 
because connectivity states arc defined by unclosed loop segments. In contrast, the final scalar product (|13p contains 
only contributions from single closed loop configurations admitting exactly one interface. 

An analog observation can be done in the plane geometry. Suppose the system is defined on a finite square of 
N = L 2 sites. A given loop configuration contributing to the partition sum Z(N\Q, r) admits one and only one interface 
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Figure 6: An example term in the matrix element T/s a of the transfer matrix. 
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Figure 7: Closing of the loop at the two boundaries of the strip. Dashed lines indicate periodic boundary conditions. 



Figure 8: Mapping of the interface problem from the plane to the cylinder. 
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Figure 9: Vertex weights in the colored loop model. 



r going from to a point of the boundary. But if one considers the subsystem contained in a smaller square centered 
on and of side L' < L, one may find several interfaces going from the origin to a point of the smaller square. 

5.4 A local model for the interface weight 

The existence of the interface T and its weight w l are non-local features, which can fortunately be expressed by local 
constraints on a colored loop model with loop fugacity n = 0. The loops may take two colors, according to the rules : 

• vertices of type (a) and (b) involve two strands of any colors, but vertices of type (c) must involve two strands 
of the same color. 

• the vertex weights are those given in figure 

• a strand must change its color when it crosses a seam going along the cylinder. 
The loop model defined by these rules satisfies the following properties : 

• Every colored loop configuration on the cylinder admits exactly one interface. 

• If a loop configuration admits a unique interface, then it can be colored in exactly two ways. 

• The interface length I is equal to the number of vertices which involve two different colors. 

Thus the partition function of the colored loop model is equal (up to a constant factor) to the partition func- 
tion Z cy i(M,L). 

5.5 Other partition sums 

In the previous discussion we described the partition sum of one loop with one interface. However, the colored loop 
model allows us to compute partition sums with up to two interfaces, and/or containing open trails. 

The trick used to impose the presence of two interfaces is simply to eliminate the color change across the infinite 
seam. 

Open trails are introduced in a standard way into the connectivity basis, by adding states where one or more points 
are connected to a "virtual" partner. The additional rule is that these points cannot be connected with each other. 
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Figure 10: Estimates of the exponent -X"(o.i) for 



Denote n\ the number of open trails and n% the number of closed loops. The cylinder circumference L must have 
the same parity as n\. The partition sums which we are able to reach with the colored loop model have n\ + ni < 2 
interfaces. The corresponding physical exponent is denoted X( ni ri2 y 

5.6 Numerical results 

Our numerical procedure is to diagonalize the transfer matrix of the colored loop model, using the power method. 
Define the exponent -X"(o,i) by : 

Z p i(N\0,r) 1 



Using equation (jTTj) , we define the finite-size estimates X( ^(L) by : 

L 2 



^(0,l)(i) = ^(/l 0) -/L,(0,l)) 



(15) 



(16) 



These estimates are also useful to determine the position of the critical point w* . Indeed, consecutive estimates 
X(n t u(L) and X( 0i i)(£ + 2) coincide at w = w*(L,L + 2), and — by analogy with standard phenomenological renor- 
malisation group methods — we expect this value to converge to w* when L is large. This method is illustrated in 
figure [TOl The resulting critical line is plotted in figure ITT1 

Following a similar method, we define the general exponent Xi ni ,n 2 ) by comparing the free energy fL,(ni,n 2 ) with 
the free energy associated with the greatest eigenvalue for the same system width : 



X 



L 1 



(n 1 ,n 2 )(I J ) — 2^ ^yf L fh,(ni,n 2 



(17) 
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Figure 11: Position of the critical point w* as a function of x. The values were obtained for width L = 8, 10. The 
curve is very well approximated by the formula : w* — 1 + 0.4 x. 

Note that ff ] is the free energy of one unconstrained loop (resp. trail) when L is even (resp. odd). 

The errorbars are given by the spacing of the consecutive intersections of the finite-size estimates. The results are 
compared with the empirical formula : 

m 2 — 1 3 

X (nun 2 ) = — ^2 — ' m = 2 Ul + 2U2 ~ 1 ^ 18 ^ 

to be discussed in the last section. 

One must consider these results cautiously, because in the loop model, logarithmic corrections appear in quantities 
such as . 



Name Description Numerical Formula (fTS)) 



X (0 ,i) 1 loop, 1 interface 0.04 ± 0.005 

X {1 [ 0) 1 trail, 1 interface 0.00 ± 0.005 -0.06 

X (0 ' 2 ) 2 loops, 2 interfaces 0.80 ± 0.02 0.67 

X(i'i) 1 loop, 1 trail, 2 interfaces 0.50 ± 0.008 0.44 

X( 2 ,o) 2 trails, 2 interfaces 0.21 ± 0.02 0.25 

Table 1: Estimates of the exponents X( ni1l2 ) by transfer-matrix diagonalisation. 
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6 Monte-Carlo simulation on a square 



In order to access directly some geometrical features of the model, we simulate it by a Monte-Carlo procedure [18] on 
an L x L square lattice. First, we describe a Metropolis procedure used for the loop model defined by the partition 
sum Zo(N). Then, we modify slightly this procedure to simulate the ensemble defined by Z(N\0,r). We discuss the 
main aspects of the algorithm : detailed balance, ergodicity and autocorrelation times. At the end, we give some 
numerical results obtained by this method. 

6.1 Metropolis procedure for the loop model 

The statistical ensemble of interest is the set of loop configurations consisting of exactly one loop on an L x L square 
lattice. Each loop configuration C is given the Boltzmann weight II(£) = x Nc , where N c is the number of crossings of 
the loop with itself. The boundary of the square consists of "frozen" vertices, as shown in figure 03 

The idea of the algorithm is to modify locally the system without violating the one-loop constraint. Let us denote 
Pij the 2x2 "plaquettc" consisting of the vertices {(i, j), (i + 1, j), + 1), (i + + 1)}, and the state of these 
four vertices. Every plaquette has eight "external legs". Outside Py, the loop connects these legs with each other. 
Fixing this external connectivity — call it a — defines the set B a of plaquette configurations (3 respecting the one-loop 
constraint when adjoined with a. 

An elementary iteration consists of four steps : 

1. Pick uniformly a plaquette P among (L — 1) x (L — 1) possibilities. The state of the plaquette is denoted (3. 

2. Visit the loop to determine the external connectivity a of the legs of P. 

3. Pick uniformly a configuration (3 1 from B a . 

4. Perform the local change (3 — > (3 1 on P with probability : 

wp^n.m '"W<n(ffl (19) 

I 1 otherwise 

This algorithm satisfies the detailed balance condition for the probability distribution H(C) on the set of one-loop 
configurations : 

V£, £ Il{£)W(£ -> £) = IL(£)W(£ -» £) (20) 

6.2 Metropolis procedure for the loop model with an interface 

Consider next the statistical ensemble of one-loop configurations admitting an interface T that connects the center O 
of the square to a point r on the boundary of the square (see figure 02) ■ Each loop configuration is given the weight 
II(£) = x c w , where N c is the number of crossings of the loop with itself, and I is the length of the interface T. The 
above Metropolis procedure can be adapted to simulate this ensemble, and is defined by the elementary iteration : 

1. Pick uniformly a plaquette P among (L — 1) x (L — 1) possibilities. The state of the plaquette is denoted j3. 

2. Visit the loop to determine the external connectivity a of the legs of P. 

3. Pick uniformly a configuration (3' from B a . This defines a new loop configuration. 
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Figure 12: Two examples of plaquette flips which alter locally the interface T (T is represented as a thin line). 



4. Test whether the new loop configuration still admits an interface. If it is the case, denote V the length of the 
new interface. If not, reject the step and skip to a new iteration. 



5. Perform the local change (3 ~ > (3 1 on P with probability 



W(/3 ->/?') = < x"o(ff) w i 11 x w x w 
] 1 otherwise 



(21) 



This algorithm is able to modify locally the interface, when the plaquette picked at step 1 is close to the interface. 
Examples of this process are given in figure 1121 



6.3 Characteristic times 

The Metropolis algorithm described above defines a Markov chain on the set of one-loop configurations with an 
interface : (£o, £i, ■ ■ ■ £/<")■ We are interested in the calculation of an observable : 



(A) = 



E £ n(£)A(£) 
E £ n(£) 



(22) 



where the sums are over one-loop configurations with an interface. The Markov chain is used to build an estimator 
for this quantity : 

1 Q_1 

^4fco,Q,T = q ^ A(C ka+qT ) (23) 
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If the configurations Ck were independent and obeyed the equilibrium distribution II(£), then the first two moments 
of the estimator Ak 0> Q >T would be : 

(A ko , Q , T ) = {A) (24) 
([A k0 ^ T -(A)] 2 ) = l([A-(A)} 2 ) (25) 

where the brackets denote averaging with respect to the distribution II(£). The parameters introduced in equation (|23[) 
play the following roles : 

• The initial time fco is the number of configurations to discard at the beginning of the Markov chain. Indeed, the 
first configurations do not obey the equilibrium distribution II(£). The value of kg must be greater than the 
equilibration time of the observable A under the process. 

• The time interval r allows us to discard redundant information arising from non-independent configurations Ck , 
Ck+i, etc |19j . One has to choose a value r greater than the auto- correlation time of the observable A under the 
process. 

• The number of samples Q controls the accuracy of the estimator. If the interval r is large enough, the error on 
(A) is proportional to Q^ 1 / 2 . 

In order to determine the equilibration time and the auto-correlation time for the obscrvables of interest, we make 
a specific hypothesis. The obscrvables we want to compute depend only on the interface. We define the "interface 
overlap function" : 

Sri cn = 2 x < rnr '> ,20) 

The hypothesis is that characteristic times for observables depending only on the interface are bounded by a typical 
timescale associated with the decay of Sr(£o, £ T ) as a function of r. This decay is well described by the auto-correlation 
function : 

Sr{C ,£ T ) - S r {C ,C K ) „ , or . 
<M r ) = , r , K->oc (27) 

This function is plotted in figure IT5T The time unit used for plotting is one Monte-Carlo sweep (mcs), corresponding 
to L 2 elementary iterations of the Metropolis algorithm. 



6.4 Ergodicity 

The Markov chain (£o, £i, ■ ■ ■ ) defined by the Metropolis procedure converges to the equilibrium distribution II(£) 
under two conditions : detailed balance and ergodicity. A Markov process is said ergodic if and only if for any pair of 
configurations (£,£'), there exists a path of transitions, going from £ to £', with non-zero transition probability : 

£^£!^£ 2 ^...£ fc ^£' (28) 

For the two versions of the algorithm, we do not have any proof of ergodicity for a system of arbitrary size. Nevertheless, 
we have been able to check exactly the ergodicity property for systems of size L = 2, 4. The transfer-matrix algorithm 
is used to count the number of allowed configurations for small width L. The Metropolis algorithm is modified to visit 
every configuration obtained by a sequence of local changes, starting from an arbitrary configuration. We find that 
every allowed configuration is reached by the process. Table [2] gives the number of one- loop configurations on a L x L 
square, for the two cases. 
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L M loop A/"i loop, i interface crgodicity checked 



2 40 32 x 

4 5 373 952 3 072 000 x 

6 4 380 037 227 741 184 1 705 236 234 240 000 



Table 2: Ergodicity check of the Metropolis algorithms for small system sizes. 
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6.5 Numerical results 

Using Monte-Carlo simulations for various system sizes L = 16, . . . , 128, with the parameters 



fc = 100 mcs 

r = 100 mcs 

1000 samples for L < 128 
100 samples for L = 128 



Q 



we estimate the average length (I) of the interface V, and its fluctuations. We obtain an evidence (see figures IT4l and 
IT5]) for the scaling laws : 

(/) = L 1/u Hi ([w-w'iL^L 1 ^ (29) 
(Al 2 ) = {l 2 )-{l) 2 = L 2 / v H 2 ([w-w*{L)]L 1 ' v ^ (30) 

where w*(L) is the "finite-size critical point", and Hi,H 2 are scaling functions. The function H\{u) is increasing on 
the real line, and the function H 2 [u) has a maximum at u = 0. We define w* (L) as the value of the parameter w which 
maximizes the fluctuation amplitude (Al 2 ). The maximum of the quantity (AZ 2 ) is expected to follow the scaling law : 

(AZ 2 ) max cx (31) 

This scaling law is confirmed by the numerical simulation, as shown in figure [TBI The numerical fit leads to : 

l/v= 1.28 ±0.018 (32) 



6.6 Scaling laws for the interface 

We interpret the scaling laws for the length I of the interface, using the scaling hypothesis for the partition function, 
close to the critical point w — w*. Suppose that the partition function Z[N\0, r) is a two-point correlation function 
in some critical field theory. In this theory, the correlation length £ is related to the deviation from the critical point 
by: 

Z(x\w-w*\- v (33) 
The function Z(iV|0,r) defines the anomalous dimension rj : 

Z ( N \°^-^^9(r/0 (34) 

where d = 2 is the space dimension, and g is a scaling function. The partition function Z(N\0, r) is also the generating 
function for the length I of the interface going from to r : 

ologw 

where the scaling functions hi, h 2 can be expressed in terms of the function g. So, in particular, the fractal dimension 
of the interface at the critical point is<i/ = l/^=1.28. 
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Figure 14: Mean length (/) of the interface, as a function of w, for x = 1. The parameters {u>*(L)} and v arc 
determined by the maximum of (A/ 2 ) (see figure ITB)) . 



T 




(to - to*) X L 1 /' 



Figure 15: Fluctuations of the length / of the interface, as a function of w, for x = 1. The parameters {w*(L)} and v 
are determined by the maximum of (A^ 2 ) (see figure [T6]l . 
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10000 - 



1000 - 




Figure 16: Maximum value of the fluctuations of the length I of the interface, as a function of the system size L. The 
linear fit gives (AZ 2 ) max oc L 2 / u , with 1/u = 1.28. 



This result can be compared with the transfer-matrix calculations of the exponents X( niin2 ). As for the case of 
polymers [TUJ [13], the thermal exponent v is related to the "two-leg exponent" X^_ 2 ) ■ 

X (0 ,2) = 2 - - (37) 
v 

Indeed, the exponent X(o,2) corresponds to the insertion of two interfaces. According to this argument, the estimated 
value of X(o,2) ( see table [1} gives a fractal dimension : 

d f = 1.2 ± 0.02 (38) 

This estimate of v is only marginally compatible with (|32p , but this was to be expected given the difficulty in estimating 
reliable error bars within the transfer matrix method. 



7 Conclusion 

In Brownian motion, the fractal dimension of the boundary is known to be Dr = 1/^ = 4/3, a celebrated conjecture 
of Mandelbrot's, established rigorously by Lawler, Schramm and Werner At the critical value of w where the 

interface is allowed to grow in our fully packed trail model, the result (|32[) is quite close to this value. Since by 
construction T is self avoiding, and since few universality classes are known for self-avoiding walks (even though the 
background of trails induces non local interactions in T) it is tempting to speculate that our v = 3/4 indeed. 

2 We note that in this case there is no microscopic duality between the non-intersecting Brownian walk and any sort of "interface", since 
the Brownian walk is not forced to fill up space. It would be interesting to see if any such interface can be defined. 
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In the context of Brownian motion, the problem closest to ours is the problem of non-intersections of packets of 
walks. We shall restrict to situations with packets made of one or two walks only. Suppose we have n\ walks and n 2 
pairs of walks, and suppose we demand that no intersection occurs between either sets, while for each pair of walks in 
the set ri2, intersections are allowed. The exponent is then known to be [5] : 

Tfl 2 — 1 

x = 2((n u n 2 ) = — ^— (39) 

where : 

m = 2ni + 3n 2 (40) 

One recovers the usual case of non-intersection exponents when n 2 — 0, n\ = L in Duplanticr's notations [3] (not to 
be confused with L the system size in this paper). 

Assuming (without clear justification) that similar properties might hold for the non-intersection exponents of 
trails, we found that the following empirical formula : 

m 2 -1 3 
x = — — — , m=-ni+2n 2 -l (41) 

gives a decent fit to our data, as mentioned in section f5. 61 

Recall that the numerical data used to obtain this conjecture are affected by very bad convergence properties. An 
analytical approach and further numerical confirmation are needed to make progress on the problem. For example, 
one may think of attacking the problem using quantum gravity methods. 

An important point concerns the value of n (the loop fugacity) that is best suited for comparison with Brownian 
motion. We chose in this paper n = mostly by analogy with the usual SAW case, but notice that formally the 
corresponding central charge is c = — 1 [7]. This might not be the best choice, since non-intersection exponents such 
as ([39]) appear in a CFT with c = [3] . On the other hand, the value c = — 1 is really a property of the "bulk" of 
dense trails, and it is not excluded that their frontiers might be described by a different CFT. In fact, it might well be 
that non-intersection exponents are independent on the value of n — after all, the dense trail models with n < 2 have 
properties which, in many respects, are independent of n (such as the fuscau exponents, which are all vanishing [7]). 
Some thoughts on the proper generalization of the lattice model needed to deal with the case n ^ are presented in 
the Appendix below. We hope to get back to this question soon. 
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Appendix : Ideas for generalizations to n^O 

In the planar geometry, when n is non-zero, several loops are allowed to cover the lattice, so the interface defined in 
section [3] is no longer unique. 

In the case n = 0, the coloring of loops with two colors was used in section [5] to express with local weights the 
existence of an interface. When n is non-zero, this trick can inspire some generalizations of the problem. These are 
best defined on a cylinder. 

Let Z co \ be the partition sum of the model defined by the vertices weights of figure [9j and the non-local weight 
n for each closed loop. Given a particular configuration C consisting of several uncolored loops, the model allows 
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several colorings of L. Intersections among the loops of C define a partition into packets of mutually intersecting 
loops. Packets may be colored independently in two ways, but loops belonging to the same packet must have the same 
color. Thus, the number of such objects appears in the partition sum : 

Z col (M, L) = J2(li w *\ ™ #1 ° ops 2 #packcts (42) 

where the sum is over uncolored loop configurations, and the first factor stands for the local weights of figure [HI 

If a change of colors is imposed when the boundary is crossed, then no packet is allowed to wrap around the 

cylinder. This is equivalent to the existence of an interface as defined in figure [8j So, like in the case n = 0, this 

model is able to express locally the existence of an interface. Unfortunately, we have not been able to define a unique 

interface, so we cannot control its length I like in the previous case. 

A special case is when x = 0. Then a packet is identical to a loop, and the partition sum becomes : 

Z col (M, L) = (2") #loops (x = 0) (43) 

Note that this model contains the Potts model as a particular case, when w = 1. 



The form (|4"2")l of the partition sum of the colored loop model suggests another generalization of the previous 
models, obtained by introducing a second non-local weight m for each packet of closed loops : 

Y(x, n,m)=J2 x* cmssinss n * loops TO #packcts (44) 
c 

where the sum is over uncolored loop configurations. 

NB : In section [31 we saw that when (n,m) = (0,2), the model with an interface is not critical. However, it may be 
possible to define a limit n — ► for the partition sum (|4"4"]) . which defines a critical model with an interface. 
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